function [verts1 verts2 solid] = GenBaseVerts

d = [0 1];
[x,y,z] = meshgrid(3*d,2*d-1,d);  % A cube
x = [x(:);2];
y = [y(:);0];
z = [z(:);0.5];

verts1 = [x';y';z'];

% [x,y,z] are corners of a cube plus the center.
X = [x(:) y(:) z(:)];
Tes = delaunayn(X)
%tetramesh(Tes,X,ones(size(Tes,1),1));


solid.Tes1 = Tes;
solid.X1 = X;

x = x + 3;
verts2 = [x';y';z'];
X = [x(:) y(:) z(:)];
Tes = delaunayn(X)
solid.Tes2 = Tes;
solid.X2 = X;

breakhere = 1;

% pt1 = [0;1;0];
% pt2 = [0;-1;1];
% pt3 = [0;1;1];
% pt4 = [0 ;-1;0];
% 
% s = [3;0;0];
% 
% verts1 = [pt1 pt2 pt3 pt4 pt1+s pt2+s pt3+s pt4+s];
% verts2 = [pt1+s pt2+s pt3+s pt4+s pt1+2*s pt2+2*s pt3+2*s pt4+2*s];
% 
% tri1.verts = verts1';
% tri2.verts = verts2';
% tri1.x = verts1(1,:); tri1.y = verts1(2,:); tri1.z = verts1(3,:);
% tri2.x = verts2(1,:); tri2.y = verts2(2,:); tri2.z = verts2(3,:);
% 
% quadcells = [1 2 3 4; 5 6 7 8; 1 3 5 7; 2 4 6 8; 1 4 5 8; 2 3 6 7];
% 
% triidx = [];
% 
% for k = 1:size(quadcells, 1)
%    t1 = quadcells(k,1:3);
%    t2 = quadcells(k,2:4);
%    triidx = [triidx; t1; t2];
%     
% end
% 
% tri1.faces = triidx;
% tri2.faces = triidx;
